A supercritical superfluid and vortex unbinding following a quantum quench 
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We study the dynamics of the relative phase of a bilayer of two-dimensional superfluids after the 
two superfluids have been decoupled, using truncated Wigner approximation. On short time scales 
the relative phase shows "light cone" like thermalization and creates a metastable superfluid state, 
which can be supercritical. On longer time scales this state relaxes to a disordered state due to dy- 
namical vortex unbinding. This scenario of dynamically suppressed vortex proliferation constitutes 
a reverse-Kibble-Zurek effect. We observe dynamics of creation of vortex anti- vortex pairs and their 
consequent motion. Our predictions can be directly measured in interference experiments 

PACS numbers: 



The understanding of order is one of the main objec- 
tives of many-body theory. Apart from topological order 
and exotic order [2], most phases can be characterized 
through the long-range behavior of their correlation func- 
tions. For example, a system of a single species of bosons 
in two spatial dimensions at finite temperature shows two 
different regimes: at low temperatures the correlation 
function of the boson operator b{x) decays as a power- 
law G(x) = {b{x)^b{0)) ~ with a scaling expo- 
nent r, at higher temperatures it decays exponentially 
(&(a;)t6(0)) - exp(-|x|/xo), with some coherence length 
xq. The transition between these two regimes, occurring 
at the critical exponent Tc = 1, is the famous Berezinsky- 
Kosterlitz-Thouless (BKT) transition Q . It has been re- 
cently observed in ultra-cold atom systems P, 01 ■ 

The power-law scaling of the quasi-superfluid state is 
due to thermally excited phonons. The gapless spectrum 
and the low dimension of the system lead to a fluctuating 
phase at all finite temperatures. If these were the only 
excitations of importance, the system would show power- 
law scaling at any temperature with the exponent r being 
proportional to the temperature T (see e.g. @). How- 
ever, vortices can drive the system to a disordered state, 
in which the correlation function decays exponentially. 
In the quasi-superfluid phase vortex-antivortex (V-AV) 
pairs lead only to a small renormalization of the scaling 
exponent. These V-AV pairs are 'confined', i.e., they do 
not separate spatially significantly. However, above the 
transition temperature, due to entropic effects and due 
to screening, these excitations become 'deconfined', i.e., 
vortices and antivortices can separate spatially. 

In this paper we study the dynamics of such a system, 
following a fast ramp across the transition, coming from 
the ordered side. We imagine two 2D superfluids that 
are coupled by a tunneling energy J±, which leads to 
a phase-locked superfluid (SF) state, in which the fluc- 
tuations of the relative phase are suppressed, while the 
total phase correlations show algebraic scaling The 
critical temperature of this state typically lies well above 
the BKT temperature Tc of the uncoupled system. We 
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FIG. 1: Illustration of the Kibble-Zurek mechanism, which 
describes ramping across a phase transition from the disor- 
dered phase, and its counterpart, the reverse-Kibble-Zurek 
mechanism. The defining feature of the latter is the dynami- 
cal suppression of vortex unbinding, which happens long after 
phonons reach steady quasi-equilibrium state. 



ask the question how the system transitions to its new 
equilibrium state after the two SFs have been decoupled 
by turniiig off J±. (The case of ID SFs has been studied 
in Refs. 0,[1, 0, [l^-) We find that on short time scales 
vortices are not important and the system shows 'light- 
cone'-like behavior in accord with prediction of Ref. ^ij 
for ID systems. Qualitatively this behavior comes from 
the fact that at distances larger than the product of a 
characteristic velocity and the time the correlation func- 
tions are not causally connected and thus decay in time 
but do not depend on the distance. Conversely at dis- 
tances smaller than this product the correlation functions 
freeze in time and depend only on distance between the 
two points. 

The state that emerges after this vortex free evolution 
is a metastable SF state. The correlation function shows 
algebraic scaling with some exponent that can be related 
to the initial temperature and the coupling energy of the 
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initial state. However, the exponent of that state can be 
supercritical, that is, the correlation function can fall off 
faster than jxl"^/"*. This state can not exist in equilib- 
rium and can be thought of as a 'superheated' superfluid. 

On longer time scales this metastable state will relax 
to thermal equilibrium phase. If the effective tempera- 
ture of the metastable superfluid is small then vortices 
do not unbind and thus the thermal equilibrium remains 
superfluid. If the metastable state is overheated then 
the relaxation to equilibrium is accomplished by vortex 
unbinding and the system becomes thermal Bose gas. 

In Fig. [T]we illustrate this process, and contrast it to 
the Kibble-Zurek (KZ) mechanism [li, [H. The latter 
refers to the case of ramping across a phase transition, 
from the disordered side. For a 2D superfluid that would 
be the thermal Bose gas phase, in which vortices are de- 
confined. If the system is then ramped across the transi- 
tion with a fast quench, some of these vortices can survive 
on a very long time scale, their recombination with anti- 
vortices is suppressed (upper part of Fig. [1]). The case 
that we consider here is ramping across the transition 
from the ordered side. The coupling Jj_ between the two 
layers suppresses both phonons and vortices. When the 
system is decoupled phonons propagate very fast creat- 
ing a metastable quasi-equilibrium state (second panel of 
the lower part of Fig. [1]). Only much later equilibration 
of vortices occur, leading to thermal equilibrium. This 
clear separation of time scales allows us to refer to this 
process as the reverse-Kibble-Zurek (rKZ) mechanism. 

We emphasize that despite the BKT transition is 
driven by thermal fluctuations, the mechanism of vortex 
or phonon creation in the process we consider can be trig- 
gered by initial quantum fluctuations. Indeed when the 
superfluids are strongly coupled together and the initial 
temperature is small the density (the momentum con- 
jugate to the phase) has large fluctuations because of 
the zero point motion. The heating mechanism of this 
system can be thought of as enhancement of this zero 
point motion following the quench. This enhancement 
of fluctuations can not be captured within conventional 
mean-field approaches. 

To model the relative phase dynamics we consider an 
XY-model \§\ and include a hopping term 
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where is an overall (Josephson) energy scale, k de- 
scribes the ratio of kinetic and potential energies. We 
can formally replace these parameters by ^qk/tt = 2Jn, 
ttQ.q/k = U (so that Oq = a/2 JnJ7, k = 'K^j2Jn/U) 
and V{t) = 2J±{t)n/ilo, which gives a coarse-grained 
representation of the Hubbard model, where the Bose 
operators have been written in a phase-density represen- 



tation. In the Hubbard model J is the in-plane hopping, 
U is the interaction energy, n is the average number of 
bosons per site, and J± is an inter-layer hopping. In the 
continuum limit the lattice size is approximately given 
by the healing length in the system. 

It is convenient to introduce the rescaled quantities 
t — Uot/h, (j) — and n — ^f^n. In terms of these, 

the equations of motion (EOMs) are 



dt 
di 



V2. 



^sin( ^fe ^'^ )+V(0/3sin;3^.(2) 



V2 



where we defined (3 = •\/27r/K. The indices ji describe 
the four neighboring sites of site i. 

We model the relative phase using a numerical im- 
plementation of the truncated Wigner approximation 
(TWA) (see Refs. [lj| for a review): The expectation 



of any quantity at some time t > can be determined 
by sampling over a Wigner distribution at time t = 0, 
and solving the classical equations of motion from to 
t. This approximation is accurate either at short times, 
or with nearly harmonic systems [l5| . In our case initial 
"light-cone" stage of dynamics is well described within 
the quadratic Bogoliubov theory where TWA is exact. 
At longer times, when nonlinear dynamics takes over, we 
expect TWA to remain applicable because by that time 
relevant momentum modes become highly occupied jus- 
tifying the validity of the approach [lj|. The advantage 
of this method is that each run strongly resembles a sin- 
gle realization of experiment, especially for a positively 
defined initial Wigner function. Thus to a very good pre- 
cision one can mimic actual experiments. In Fig. [2] we 
plot time evolution of the relative phase for a single real- 
ization on initial conditions. One can observe V-AV pair 
formation at short times and their consequent unbinding. 

We can calculate the Wigner distribution at t = 0, as- 
suming that J± is larger than the other energy scales. In 
this limit the phase fluctuations can be described within 
the Bogoliubov approximation, so the system reduces to 
a sum of oscillators. The Fourier modes 
t = are distributed according to (see Ref 
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with a = l/y^SaJ^, = 1/ tanh((jjq/2ro), and ujq = 
y/A sm{qx/2)'^ + 4 sin(9y/2)2 + VP"^, Tq being initial tem- 
perature. We use this method to extract the equal time 
correlation function: 

G(x, t) = (exp[i V20J (t) - I V20J+, (t)]), (4) 

where x is an integer separation between the points and 
t is the time after decoupling (see Fig. [3]). Because we 
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FIG. 2: Temporal evolution of the relative phase for a single realization of initial conditions. The parameters of the systems 
are V — 100, k = 10 and To = 2. Vortices are marked red, anti-vortices blue. 




FIG. 3: Plot of short (left) and long (right) time behavior 
of the correlation function as a function of time and space. 
The parameters a To = 3, for k = 10 and V — 20 for the 
left graph and To = 1, k — 8, and y = 80 for the right 
graph. At short times the dynamics separates into instan- 
taneous, damped oscillatory behavior, and a 'light cone' like 
pulse forming metastable quasi-superfluid state. At longer 
time scales the correlation function shows exponential decay 
due to dynamical vortex unbinding. 




are using periodic boundary conditions G{x,t) depends 
only on the separation between the points x. Note that 
this correlation function (or rather dx' G{x' , t\) can be 
directly measured in interference experiments [l|, 0, [l3| ■ 
We indeed see very clear emergence of the light cone ther- 
malization: At separations larger than 2vt, where v is 
characteristic phonon velocity, G(a;, t) is almost x inde- 
pendent - it uniformly decreases in time. Once 2vt > x 
the correlations freeze in time and depend only on x. We 
find that the state that emerges within the light cone 
shows algebraic scaling, and therefore can be referred to 
as a quasi-superfluid. At much longer times scale the cor- 
relation function clearly relaxes to the exponential equi- 
librium shape due to vortex unbinding. 

Due to the relation between the power-law exponent 
and temperature in equilibrium, we can associate an ef- 
fective temperature T* with the state inside the light 
cone. We can estimate this effective temperature by con- 
sidering the linearized dynamics, which neglects vortices. 
Then Eqs. ^ have the solution: 

0fe(<) = ^fc(O) cos cjfci - rife (0)/ti;fe sin cjfci (5) 
nfe(t) = (i>k{0)uJkS\\iujkt + nk{Q)cosujkt. (6) 



We then consider the quadrature 



{<f>l{t)) = 



2uj 



k,Q 



■ cos^(wfci) -I- 



■ sin^ (ujkt). (7) 



FIG. 4: a) T* , as given in Eq. (8] for «: = 1,3, 10, from top 
to bottom, and for V = 20. The line Tc = tv/2 is the BKT 
temperature of the uncoupled SFs. b) - d) Simulations for 
these values of k and V. 



where rk^o and cuk.o correspond to the quantities at i = 0, 
that is to V{0) = V, and ujk to the dispersion with 
V — 0. The long-time limit {(j)l.{t — > oo)) is rkfi/4:UJk,o + 
rk,o^k,o/4^i^k- We equate this formally to a thermal en- 
semble of momentum-dependent 'temperature' T^, i.e. 
{(j)l{t oo)) = rl/ujk, where r* = 1/ taiih{LUk/2T*). We 
solve for T^, and find that for large it simplifies to 
a single value, independent of k: 



rj-Tk 



4tanh(yF;52/2ro)' 



(8) 



For small To we have T* 



2J±/J in terms 



of original Hubbard parameters), that is, the temper- 
ature is fully determined by the initial coupling energy. 
The coupling energy between the two layers is transferred 
into the in-plane kinetic energy. We note that if J± ex- 
ceeds the chemical potential Un the quantum rotor model 
becomes inadequate and then T* saturates at the value 
T* - Un/J. For large Tq we have T* « ^. This is a 
reflection of the doubling of the degrees of freedom when 
two layers are uncoupled. In Fig. 2] a) we plot T* for 
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FIG. 5: Exponent r extracted from fitting tiie correlation 
function G{x) to the algebraic form, for different initial cou- 
plings. We use To = 1 and k = 8, and V — 80, 70, 50, and 
20, corresponding to curves I to IV. Except for the curve I 
corresponding to Fig. |3]-b) the correlation function can be fit- 
ted with an algebraic function throughout the time interval 
shown here. 



V = 20 and for different value of k. For k = 1 cor- 
responding to large J±, T* is always above the critical 
temperature Tc = 7r/2, for k = 10, it crosses it as initial 
temperature increases. We therefore expect to see very 
little vortex formation for small initial temperatures for 
K = 10, and many vortices for all temperatures for n — 1. 
This is indeed the case as we show in Fig. 01 which plots 
the number of vortices riy versus time and initial tem- 
perature for three different values of k. This number is 
obtained by counting the plaquettes with a phase wind- 
ing of 277, and then by averaging over many runs. 

To quantify the crossover from the supercritical super- 
fluid state to the normal phase we fit the correlation func- 
tion G{x, t) at different times numerically using algebraic 
{G{x,t) = c{L/'K\sm{Trx/L)\)~'^/'^) and an exponential 
{G(x,t) — cexp(— I sin(7ra:/L)|/a)) functions. In equilib- 
rium the algebraic exponent r would be the relative tem- 
perature T /Tc- Any value above 1 is therefore supercrit- 
ical. We use these two fitting functions in four examples 
with To = 1 and k = 8, but with different initial cou- 
plings V = 80, 70, 60, 20. In Fig. [Slwe show the exponent 
r as a function of time for these cases. In all cases G{x, t) 
develops algebraic scaling after the light-cone dynamics. 
Note that for large initial couplings (I-IH) the emerging 
scaling exponent r is well above the critical exponent. At 
longer times the exponent slowly increases in time. Dur- 
ing this process, the correlation function can still be well 
fitted with an algebraic function. Eventually the correla- 
tion function develops exponential scaling. This regime 
is reached for I^ = 80 (I) within the time interval shown 
in Fig. [HI signalling that the thermal Bose gas phase has 
been reached due to vortex unbinding. For = 70 (II) 
and V = (III) the time scale of vortex unbinding is 
longer then the time interval shown. For V = 2Q (IV) the 
system remains in a quasi-superfluid state. We conclude 
from these examples that the supercritical superfluid has 
a long life-time and the exponent r can significantly ex- 



ceed the maximum equilibrium value. Such metastable 
states should be thus experimentally feasible. 

In conclusion, we have studied the dynamics of the rel- 
ative phase of a bilayer of superfluids in 2D, after the hop- 
ping between them has been turned off rapidly. We find 
that on short time scales the dynamics of the correlation 
function shows a "light-cone" -like behavior. The system 
then develops an algebraic phase that can be stable or 
metastable. The latter can be thought of as a super- 
heated superfluid. On long time scales this metastable 
superfluid relaxes to a disordered state via creating of 
vortex-antivortex pairs and their consequent unbinding. 

Our predictions can be directly probed in experiment. 
The behavior of the relative of phase of two superfluids 
can be studied by interference experiments 17|, so both 
the supercritical state and the algebraic-to-exponential 
evolution can be tested. Vortex unbinding can also be 
demonstrated by direct observation of free vortices jlj. 
The time scale of the light-cone dynamics can be esti- 
mated as tLC ^ L/v^ where the system size L ^ 10~'*m 
and the phonon velocity v ~ 10~^m/s leads to ^lc ^ 
0.1s. The rate of vortex unbinding is exponentially sup- 
pressed compared to the BKT energy scale [l8| . 
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